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METHOD AND SYSTEM FOR ADAPTIVE 
CHANNEL ESTIMATION TECHNIQUES 

FIELD OF THE INVENTION 

In general, the invention relates to wireless connmunication systems. 
More specifically, the invention relates to cellular communications and in 
particular, to a method and system for channel estimation within cellular 
communications. 

BACKGROUND OF THE INVENTION 

For optimal reception of a digital signal over a broadband cellular 
channel, the channel to a desired transmitter must be determined. Because 
channels can consist of a number of multipath components or rays each with 
its own time of arrival and Doppler frequency, a channel estimator should 
determine the actual time of arrivals (TOA's), Doppler frequencies, and their 
corresponding complex gains. The channel estimator should also adapt to 
many different transmitting environments where there are a very large number 
of rays with a TOA spread and a Doppler frequency spread. A maximum 
likelihood (ML) channel estimation approach as is known in the art. may 
provide the best TOA and Doppler frequency estimates, however it lacks a 
closed form solution. 

Prior art TOA estimators such as a TOA estimation using ESPRIT, are 
not robust to channel conditions with a large number of rays with a TOA 
spread and are not suitable for multi-user channel estimation. Other channel 
estimators may assume a worst-case channel and thus need to estimate 
many more parameters than is necessary under typical channel conditions. 

Thus, there is a significant need for a method and device for improving 
multi-user channel estimation techniques that overcome the above 
disadvantages and shortcomings, as well as other disadvantages. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an overview diagram of a preferred embodiment for a cellular 
communication system In accordance with the present invention; 
5 FIG. 2 is a block diagram illustration of a communication device for the 

cellular communication system in FIG. 1, in accordance with the present 
invention; 

FIG. 3 is a block diagram illustrating details of a channel estimation 
device for the communication device of FIG. 2, in accordance with the present 
10 invention; 

FIG- 4 is a flow chart representation of a method performed by one 
embodiment of the Channel Estimation Device of FIG. 3. in accordance with 
the present invention; 

FIG. 5 is a flow chart representation of a method performed by a 
H 15 second embodiment of the Channel Estimation Device of FIG. 3. in 

nj 

M accordance with the present invention; 

% FIG. 6 is a flow chart representation of a method performed by a third 

K embodiment of the Channel Estimation Device of FIG. 3, in accordance with 

the present invention; and 
20 FIG. 7 is a flow chart representation of a method performed by a fourth 

embodiment of the Channel Estimation Device of FIG. 3, in accordance with 

the present invention. 
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DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 

FIG. 1 illustrates a wireless communication system 100 in accordance 
25 with one embodiment of the present invention. As shown in FIG. 1, a base 
station 110 provides communication service to a geographic region known as 
a cell 103. At least one user device 120 and 130 communicate with the base 
station 110. 
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As shown in FIG. 1, user devices 120 have a single antenna 101, while 
user devices 130 have at least one antenna 101. One embodiment of the 
invention provides that the user devices 120 and 130, as well as the base 
station 110 may transmit, receive, or both from the at least one antenna 101. 
An example of this would be a typical cellular telephone. Additionally, one 
embodiment of the invention may provide that the communication system 100 
uses techniques such as orthogonal frequency division multiplexing (OFDM). 
Furthermore, one embodiment may provide that alternative transmitters and 


Q 10 receivers known in the art may be used and user devices as well as base 

O 

M stations may be interchanged and referred to commonly as communication 


devices, transmitting units, receiving units, transmitters, receivers, 
UJ transceivers, cell phones, or any like term known in the art. 

g Illustrated in FIG, 2 is a block diagram for one embodiment of a 

^ 15 communication device 200 that may be used within the wireless 
H communication system. The communication device 200 may include at least 

p one antenna 101 wherein the outputs of the antenna may be provided to a 

^ receiving unit 201. The outputs of the receiving unit 201 may be provided to 

at least one antenna combiner 202. The antenna combiner 202 may process 
20 multiple signals into one signal. The signals from the receiving units 201 may 
also be fed into a combiner controller 210, which can regulate the operation of 
the at least one antenna combiner 202. The signals from the receiving units 
201 may also be fed into a channel estimation device 208. A pilot symbol 
generator 212 capable of generating pilot symbol information may be used by 
25 the combiner controller 210 to control the antenna combiner 202. The pilot 
symbol information generated by the pilot symbol generator 212 can also be 
used by the channel estimation device 208 to estimate the channel responses 
between a transmitting device and the communication device 200. The output 
of the antenna combiner 202 may be fed into an information decoding unit 
30 206, which can decode an antenna combiner output 204 and generate data 
information 21 3 that was received by the antennas 101. 
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A channel estimation device 300 suitable for use in one ennbodiment of 
the communication device is illustrated within FIG. 3. One embodiment of the 
channel estimation device 300 may include a received training memory unit 
5 310, capable of storing output signals from a receiving unit. A microprocessor 
340 may collect the pilot sequence information from a pilot sequence input 
350 and may receive inputs from a TOA matrices memory unit 320. The 
received data stored in the received training memory unit 310, and the inputs 
from the microprocessor 340 may be fed into a matrix multiplier 330 where 

10 they are multiplied. The result can be used in iteration calculations to 

compute a channel estimate in one embodiment of the invention. In another 
embodiment of the invention, the TOA matrices memory unit 320, the matrix 
multiplier 330, and the microprocessor 340 may be implemented in a 
processor 380. Additionally, an embodiment of the invention may provide the 

15 microprocessor 340 to make use of a Fast Fourier Transform device 390 to 
efficiently implement the channel estimation algorithms of the invention. 

FIG. 4 is an illustration of a block diagram for a frequency-domain 
multi-user channel estimation method 400 in accordance with the invention. 
In order to describe the method of frequency-domain multi-user channel 

20 estimation 400, multi-user channel estimation based on time of arrivals 
(TOAs) for multi-path components of a transmitted signal will first be 
described. 

Using prior art TOA estimators such as estimation of signal parameters 
via rotational invariance techniques (ESPRIT), TOA's for use in OFDM 

25 channel estimation may be found; however, ESPRIT fails when there are 
many multipath arrival times, or when the communication system is multi- 
user. In addition, ESPRIT is not an optimal TOA estimator since it is not a 
maximum likelihood (ML) estimator. The performance of ESPRIT may also 
suffer when the pilot symbols are not constant modulus. This can occur when 

30 the received signals need to be divided by the pilot symbols, which can result 
in an emphasis on subcarriers where the pilot is weaker in power. 
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In one embodiment of the invention, a near-ML technique for multiple 
sinusoidal estimation may be extended to a multiple user TOA/channel 
estimation in OFDM or for any frequency-domain channel. The received 
OFDM signal in the frequency-domain may be given as: 


(Mxl) Y(A:) = Xh„(^KW + N(^) (1) 

where M is the number of receive antennas, U is the number of users or 
transmitters, Hu{k) (an Mx^ vector) is user u's channel on subcarrier /c, Xu{k) is 
user i/'s symbol on subcarrier /c, and N(/c) is a vector of additive noise that 
may have covariance matrix a^l where cr,^ is the noise power. The 

dependence on time (or baud number) has been removed to simplify the 
detailed description. User u's channel on antenna m may be modeled by the 
equation: 

^.,«W = f/'(^)y™,„ (2) 
where the P^xl vector, y^ ^ , is user u's channel gains on antenna m, Pu is the 

number of distinct TOA estimates for user u, and the Pyx1 vector fu{k) may be 
Represented by the equation: 

(3) 


In equation (3), tu,e is the scaled TOA (i.e., not the actual TOA) for user u and 
path e (^=1,...,P(/). The true TOA's can be recovered by multiplying f^,^ by 
K/{2nW) where W is the bandwidth of the K subcarriers. 

For one embodiment of the invention, equation (2) is essentially a sum 
of weighted sinusoids. The sinusoidal frequencies of equation (2) (which 
correspond to the TOA's) and their relative gains need to be determined to 
establish the channel estimation. Also, the TOA's corresponding to a single 
user may be assumed to be the same across the receive antennas. If 
hardware implementations cause the synchronization to be different on each 
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receive antenna or cause the delay to be different due to filtering effects, then 
the TOA estimator may need to be run separately on each receive antenna. 

Stacking the received frequency-domain signal on antenna m across K 
subcarriers may be expressed as: 

= I;XAy«,.+N„ (4) 


jyhere KxK Xu and KxPu Fy are: 


i/=i 


f„''(0) 


(5) 


X„ =diagK(0),...,x„(/^-l)} and F„ = 

f/' (/:-!) 

Within another embodiment of the invention, the goal of channel estimation 
may be to find for each user the number of TOA's, Pu, the TOA's themselves, 
tu,e. and their gain values, y„,^ . A multi-dimensional search over all of these 

unknown parameters to find the best fit according to (4) may be a prohibitively 
complex task. Therefore the problem can be separated into two more 
tractable problems where the gains can be eliminated from the estimation 
problem and then a search can be performed just over the number of TOA's 
and the TOA's themselves. To determine the gains, a Least Squares (LS) 
approach may be taken with a weighting, a, placed on the magnitudes of y„,^ 

so that the inverse in the following equation (7) is not ill conditioned. For the 
following equation (6), when a=0 the solution is the ML estimator. The 
previously described parameters may be found as the solution to the 
equation: 


M 


min 
t P y ^ 


V — X F V 

m X ^ u u * m.u 


M U 


(6) 


The solution of (6) for the magnitude vectors, y^,^ , given the TOA estimates 
can be shown to be: 


' mJU 
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where P = ^P^ , and PxP R and Px1 are given as: 


t/»i 


R = 


R 


and H„. = 


F-XrY„, 


(8) 


and Pu^Pv Ru,v is given as: 

R„,, = F„^Xf X,F,. +^^(1/-. v)!,,,^^ (9) 

where l^xm is an nxm matrix of all zeros except the upper left 
min(/7,m)xmin(/7,m) corner which is an identity matrix of size 
min(n,/77)xmin(n,n7) and 5(/?)=1 when n=0 and 5(n)=0 when n^O. 

The mean squared error of the channel estimator can be shown to be: 


(10) 


10 The error does not depend on y„,^ and thus the TOA's can be searched 

independently of the gain values. This leads to the near-ML TOA search 
described below. At each step in the algorithm, a two-dimensional search 
may be performed to find an additional TOA and the user to which it belongs 
so that the MSE in equation (10) can be minimized. Because only a single 

15 TOA is added during an iteration of equation (10), R'^ may be updated at each 
step by using the rank one inverse update formula (the rank one inverse 
update formula is known in the art). By only adding a single user's TOA at 
each step, the algorithm can accommodate the case of each user having a 
different number of TOA's. 

20 The following is the near-ML multi-user TOA search, and [] is used to 

represent an empty matrix. First, let Kxl Gu,m and Kxl Zu{t) be: 


G„,„=xrY„ and z„(0 = 


1 


(11) 
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1. Eo=TK\' . ^=0. Ro' =□. ro=n. Fo=[]. Ho=[]. qo=D. and 

i»i=l . - 

du=cc + ^\x„ ikf for i/=1 , . . . ,U. 

M 


i=0 


2. £=e*A, (^-l)xlr, =F,^,X.,2„,(/,). E,=E,_,- 


r 

where (^xl) q,„, = 


(-\'t 


.1 


J 


Perform a two dimensional search over t 


and {Uf is which user's TOA minimizes E^) and find the and Ug to minimize 
and then update: 


(^xl)H„,,,= 


-1 


0"^ 0 


-r;>. 
1 


, F, =[F,. J X„2„^ (/,)], and 
l-rl'Rt, l] 


3. If stop criteria as described below is not reached then go back to 


10 step 2. 

4. The total number of TOA's, P. is the TOA's are fi through f^, 
and the users associated with each TOA are u^ through u^. For one 
embodiment of the invention, the search over TOA's, f^, performed in step 2 
15 may be done over some finite set of times from fmin to fmax- Since there are a 
finite number of TOA's to search over, and G^„,z„^ (t^) can be pre- 

computed. 

In one embodiment of the Invention, the stopping criteria may be to 
stop when the error at step £, E^, is close to the noise floor. This embodiment 
20 of course, requires an estimate of the noise power. Another embodiment of 
the invention may define the stopping criteria as "stop if the value of the Peak 
to Average Power Ratio (PAPR) of the error update at step £ versus all of the 
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TOA's being searched is less than some predetermined threshold". The error 

X . H e.ni" t,m 


update at step i may be expressed as the equation = 


An example of a preferred embodiment may set the predetermined 
threshold to five. The PAPR stop criteria can provide a value for crl by 

5 dividing by MK. An example of this is the mean squared error averaged 
over all M antennas and all K subcarriers. Additionally, once the TOA*s are 
found, the channel gains, y„,, can be found using equation (7). 

One embodiment of the invention may extend the previous algorithm to 
both TOA and Doppler frequency estihiation. The received OFDM signal in 
1 0 the frequency-domain may be expressed as the equation: 

(Mxl) Y(/:,6) = XH„(^,Z>K(^,6) + N(/:,6) (12) 
if=i 

\ ; " where M is the number of receive antennas, U is the number of users or 

=P transmitters, Hu{Kb) is user u's channel on subcarrier k on baud number b, 

P 

15 XuiKb) is user u's symbol on subcarrier k and baud and H{k,b) is a vector of 
additive noise with covariance matrix all . User u's channel on antenna m 
may be modeled by the equation: 

= C(^,^^)y.,„ (13) 
where P^xl y„, „ is user u's channel gains on antenna m, Pu is the number of 

20 TOA and Doppler frequency estimates for user u, and P^xl fu{k,b) may be 
expressed as: 


(14) 


where tu,e and /w,^ in (14) are the scaled TOA and Doppler frequency 
respectively, for user u and path e (1<^<Pu). Within one embodiment of the 
25 invention, the true TOA's can be recovered by multiplying tu,e by K/{2nW) 

-9- 
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where Wis the bandwidth of the Ksubcarriers. Additionally, the true Doppler 
frequencies can be recovered by dividing Z^,^ by 2nANs where A is the time 
between OFDM time-domain samples and Ns is the number of time samples 
between two adjacent OFDM bauds such as the number of samples between 
bauds b and 

For a fixed b, equation (13) is essentially a sum of weighted sinusoids 
just like the TOA-only embodiment. For channel estimation within another 
embodiment, the sinusoidal frequencies across frequency (these correspond 
to the TOA's), the sinusoidal frequencies across time (these correspond to the 
Doppler frequencies), and their relative gains may need to be determined. 
Further, the TOA's and Doppler frequencies may be assumed to be the same 
across the receive antennas for a single user. If hardware implementations 
cause the synchronization to be different on each receive antenna, cause the 
delay to be different due to filtering effects, or cause the frequency offsets to 
be different, then the TOA/Doppler frequency estimator may need to be run 
separately on each receive antenna. 

Another embodiment of the invention may assume that there are B 
bauds containing all known pilots on all K subcarriers. The extension to 
arbitrary pilot formats is straightfonA/ard to anyone trained in the art. Stacking 
the received frequency-domain signal on antenna m across K subcarriers and 
B pilot bauds gives the following equation: 

= ZX.F„Y.,„+N„, (15) 


n,(o,o) 


(KBxl) Y„,= 


Y„iK-l,B-l) 


- 10- 
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j^yhere BKxBKXu and BKxPu Fu are: 

X„ - diag{x„ (0,0), . . . , x„ (/: ^ 1,0), x„ (0,1), . . . , x„ (/: - 1, 5 - 1) } 
C(0,0) 


and F„ = 


f/'(/:-i,o) 
C(0,i) 

f„'^(ic-i,^-i) 


(16) 


For one embodiment of the invention, the goal of channel estimation 
may be to find for each user the number of TOA's/Doppler frequencies, Py, 
the TOA's and Doppler frequencies themselves. tu,e and /y.^, and their gain 
values, y^j . To determine the gains, a Least Squares (LS) approach may be 

taken with a weighting, a, placed on the magnitudes of y^^ so that the 

inverse in the following equation (18) is not ill conditioned. It is obvious that 
when a=0 the solution to the following equation (17) is the ML estimator. The 
parameters are found as the solution to: 


mm ^ 


M U 


(17) 


The solution of equation (17) for the magnitude vectors, y^, ^ , given the TOA 
and Doppler frequency estimates can be shown to be: 


r 


y m,U 


(18) 


where P = > ^ ^1 are given as: 


R = 


\JU 


R 


u,u 


and H„ - 


(19) 


and PuxPv Ru,v is given as: 


(20) 
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The mean squared error (MSE) of the channel estimator can be shown to be: 


of the gain values. This leads to a near-ML TOA/Doppler frequency search. 
At each step in the algorithm, a three-dimensional search may be performed 
to find an additional TOA, an additional Doppler frequency, and the user to 
which they belong so that the MSE in equation (21) is minimized. Because 
only a single TOA/Doppler frequency pair is added during an iteration of 
equation (21), R"^ can be updated at each step by using the rank one inverse 
update formula, as is known in the art. By only adding a single user's 
TOA/Doppler frequency pair at each step, the algorithm can accommodate an 
embodiment of the invention where each user may have a different number of 
TOA's and Doppler frequencies. The following procedure may provide the 
near-ML multi-user TOA/Doppler frequency search, where Q is used to 
represent an empty matrix. First, let BKxl Gu,m and SKxl Zu{t,f) be: 



In one embodiment of the invention, the error may not depend on y^j 
and thus the TOA's and Doppler frequencies can be searched independently 


^-7(0/-0/) 


X^Y^ and z„ (/,/) = 


^->((/:-i)/-o/) 


(22) 


M 


1. ^0 = ZlY^r . ^=0, Ro=[], ro=[l. Fo=n, Ho=[l. qo=[]. and 


111=1 


d„=a^YT.V>.ik,bf foru=1 U. 


2. 


^ 1 



Where (^xl) q,^, = 
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Perform a three dimensional search over t^, f^,and (which user's 
TOA/Doppler frequency pair minimizes E^) and find the f^, f^. and to 
minimize and then update: 


(^xl) H„,,,= 


step 2. 


0"^ 0 


. =[Fe-, |X,z„^(r,,/,)].and 
[-rfR- 1] 


3. If stop criteria as described below is not reached then go back to 


4. The total number of TOA's and Doppler frequencies P is the 
TOA's are t^ through f^, the Doppler frequencies are f^ through f^, and the 
users associated with each TOA are u^ through u^. 

For another embodiment of step two. the search over TOA's. t^, may be 
performed over some finite set of times from tmm to fmax and the search over 
Doppler frequencies may also be performed over a finite set of frequencies 
from fmin to fmax. Since there are a finite number of TOA's and Doppler 
frequencies to search over, and G^„,z„^ (h^fe) be pre-computed. 

One embodiment of the invention may provide for the stopping criteria 
to stop when the error at step E^. is close to the noise floor. That 
embodiment of course, requires an estimate of the noise power. Another 
embodiment of the invention may provide the stopping criteria to quit when 
the Peak to Average Power Ratio (PAPR) of the error update at step £ versus 
all of the time of arrivals and Doppler frequencies being scanned is less than 
some predetermined value. The error update at iteration £ may be expressed 

~ //„ I — • Dividing by MBK, a value for may be obtained. 

-1 
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Once the TOA's and Doppler frequencies are found, the channel gains, y^, , 
can be found using equation (18). 

Common among the embodiments of channel estimation methods in 
accordance with the present invention, is that each method adapts (is 
adaptive) to the channel conditions, thereby using a smaller set of channel 
parameters than existing channel estimators. The channel conditions may 
include multi-path components, channel gains of multi-path components, time- 
taps, maximum delay spread, directions of arrivals, and/or additional channel 
attributes as are known in the art. The adaptivity of the invention may be in 
part or whole, a function of channel estimation iteration, channel estimation 
TOA, and/or least squares channel estimation. Each method of adaptivity 
described within the detailed description of the invention may provide an 
adapted channel estimate in place of the received channel estimate thereby 
improving performance in most channel conditions. One such embodiment 
refers again to FIG. 4, where the frequency-domain multi-user channel 
estimation method 400 can provide multi-user TOA and Doppler frequency 
near-ML estimations. 

The channel estimation method 400 may be provided a noisy channel 
estimate for at least one user Gu,m, and an update vector Zu(f,0. 410 and may 
be represented by the equation: 

-JiOt-O/) 


G„.„, =X„^Y„ and z„(/,/) = 


-j((A'-l)/-(fi-l)/) 


The iteration variables 


may next be initialized that may include setting the iteration number ^, to zero 
420. The iteration variables and how they are initialized may be represented 

by the equations E,=Y\'^„f , Ro=[]. ro=[l. Fo=[]. Ho=[]. qo=[l, and 

=^-*-XSK(^»*)r for u=7,...,L/, where [] indicates an empty matrix as is 

- 14 - 
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known in the art. After initialization, the iteration nunnber may be incremented 
by one 430. The channel estimation method 400 may then calculate an error 
update at the incremented iteration for all TOA values, Doppler values, and 
users being scanned as a function of the iteration variables and the noisy 
channel estimates for at least one user 440. In one embodiment of the 
invention the TOA's t may be scanned from fmin to fmax. the Doppler 
frequencies f may be scanned from ^min to fmaxi 3nd the usefs, u, may be 
scanned from one to U. The error update at iteration £ can be a function of 
TOA. Doppler frequency, and user and for this embodiment may be equated 

M 

e it f ~ ff"=l 

as ''-«-'-"«^'V,,^here(^-^>^^'-'=^'-'X-«^-« and 


Cr ^ Z (t f ^ 

- ^'-^^"^-I^Determining which t, and u that can maximize 

the error update at iteration i can be performed with the determined values 


labeled f^, f^, and respectively 450. The error at step £, E^, may be updated 
to equate E^=E^-i-e^(f^/^,u^). Next, the following iteration variables may be 


r 


updated 460: 


0^ 


1 


^ F,=[F,_,|X„z„,(/„/,)] 


If the appropriate stopping criteria have been reached 470, the channel 
estimation method 400 may be terminated 480. If the appropriate stopping 
criteria have not been reached 470, then the channel estimation method 400 
may proceed to block 430. 

In another embodiment of FIG. 4, directions of arrivals can be scanned 
in a similar manner as the TOAs and the Doppler frequencies. Anyone skilled 
in the art can modify the channel estimator of FIG. 4 to scan for direction of 
arrivals in a near-ML sense. 
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Another embodiment of the invention illustrating an adaptive channel 
estimator is FIG. 5, where a TOA with a phase shift training channel 
estimation method 500 may provide channel estimation by shifting multiple 
user channels to non-overlapped time regions. Given that a first user's pilots 
are p^{k), then a subsequent users* pilots (say the n^^) may be determined as: 


where Ls may represent the maximum length in number of time samples of all 
users' channel responses. The phase shift training channel estimation 
method 500 may begin by being provided a noisy channel estimate and a time 
separation value 510. The time separation value Ls may be at least as long 
as the expected delay in sample times. The noisy channel estimate may be 
equated as H„(A:) = Yik)/ p^{k) . The TOA's may be estimated next 520 from 

the noisy channel estimate using any existing TOA estimator known in the art. 
The use of the training in (23) may turn the multi-user TOA estimation 
problem stated previously in equation (6) to a single-user equivalent problem 
where a user u's TOA's can be from {u-^)LsAN/K to uLsAN/K, and A may be 
the time between OFDM time-domain samples. The complexity is similar to 
the multi-user embodiment of FIG. 4 because the number of TOA's that need 
to be searched at each iteration using the phase shift training may be the 
same as before. Thus, the final step of the phase shift training channel 
estimation method 500 may be to separate the TOA's for each user given 
what time band they are in 530. For one embodiment of the invention, user 
u's TOAs are from {U'1)LsAN/K to uLsAN/K. The phase shift training can 
enable the use of ESPRIT as a multi-user TOA estimator because the phase 
shift training problem can turn the multi-user TOA search into a single-user 
equivalent. 

In another embodiment of the invention, a gradient search may be 
performed to improve the TOA estimates from any existing TOA estimator. 
This embodiment may provide that for a single user, the TOA estimation can 
be expressed as follows if the number of TOA's. P. is known: 


-j2fA{n-\)LJK 


(23) 
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min ^ 


(24) 


"where the (/Wx1) noisy channel estimate Hn{k) is H„(/:) = \{k)/x{k) and is 

a /lfx1 matrix of channel gains for the TOA. A gradient of equation (24) 
with respect to can be shown to be: 


Vr,=glm{^e>*''YfH„(^))-i;im 


. (25) 


Given the PTOA's, the gains (i^x M)y = [Y,,...,y^]^ may be found as 


. ( 


y = R-^p where R = 2f(A:)f (/:), p = (A:)H^(A:), and f (A:) = 
26) 

Consequently, one embodiment of the invention may provide a 
gradient-based search for the TOA's as follows: 

1. Start with PTOA's from an existing algorithm (e.g., ESPRIT or 
the near-ML TOA algorithm), 

2. Determine y from equation (26). 

3. Calculate Vf, for ^=1,...,P from equation (25), for £=1,....P, t^tr 
aVt^ where a is the step size. 

4. If the stopping criteria are not reached, such as a maximum 
number of iterations, then go to step 2. 

FIG. 6 Illustrates another embodiment of the invention using the 
gradient-based search for TOA's. The gradient TOA channel estimation 
method 600 may be provided a number, P, TOA estimates for ^=1 ...P, and 
a step size a 610. A noisy channel estimate may be determined as the 
equation Il„(k) = Y(k)/x(k), The TOA gradient for each TOA estimate may 
next be calculated as a function of the TOA estimates and the noise channel 
estimate 620. The TOA gradient for the TOA, f^, can be shown as the 


m 


V 


A'-l 


equation^ Vr, = 2lm{yte^*''yf H^CA:)}-^!^ 
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,1 


Ty"ypke 
p=i 


-Ml. -I,) 


, where 


~i 


The 


(PxA/)Y = [Y,,...,ypf .and 

T 

Y = R-'p where R ='2^f(A:)f p = ^f(fc)H^(it), and f (yt) = 

L 

TOA's can be updated for ^=1 . . .P using f^=f^-aVf^ where a is the step size 
630. If a stopping criterion is not reached 640 such as a maximunn number of 
iterations is not reached, the iterations may continue from block 620, 
othenvise the algorithm may terminate 650. The gradient search may not be 
the preferred method to determine TOA's, but for one embodiment of the 
invention, the gradient search may refine estimates that are reasonably close. 
This is due to the objective function having many local minimums and that the 
gradient search does not provide a good procedure for determining P and the 
initial TOA's. 

As is illustrated in FIG. 7, another adaptive channel estimation 
embodiment of the invention can include an improved least squares channel 
estimation method 700. Assuming that there is only a single user, the least 
squares (LS) channel estimator can parameterize a frequency-domain 
channel as the weighted sum of L time taps, or equivalently the A/-point DFT 
of an L-tap time-domain channel as follows: 


L-\ 
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(27) 
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where N is the DFT size of the LS estimator (in OFDM, N>K). The LS 
estimator for receive antenna m can be put in a matrix form as represented by 
the equation: 

h..= =QY„ (28) 
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where Q is LxK matrix that may provide a LS time-tap channel estimate as is 
known in the art, and Ym is defined in the previous equation (4) with U=^ . 

In sparse channels the LS estimator can be improved by not estimating 
gains for taps below a certain power level, thus improving the pilot-to- 
parameter ratio. In this embodiment of the invention, criteria may determine 
which taps are not estimated. Assuming that only white noise is present, the 
noise floor of the LS estimator can be represented by the equation: 

a, =tT„Qd (29) 

where d is an Kx1 vector of all ones. 

The criteria within one embodiment of the invention for the removal of 
taps may be to remove all taps with power below 101og,o(oi^a;,) + 5dB . Once 

the taps are removed, the gains can be found using a LS procedure similar to 
equation (7). Using the described process, the least squares channel 
estimation method 700 may begin by being provided a time-domain channel 
estimate for at least one user and a threshold value 7710. In one 
embodiment of the invention, the threshold value may be represented by the 
equation T = 101ogio(a^<r^) + 5dB . where = cr„Qd , d is a Kx1 vector of all 

ones, cx^ is the noise power, and Q is a matrix that is known in the art which 
can provide a LS time-tap channel estimate. The time-domain channel 
estimate may be provided from any existing channel estimator. In the Least 
Squares (LS) embodiment, the time-domain channel estimate on receive 
antenna m may be given as;^ 
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h.. = 
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= QY„. , where 
/i„(I-l)J ^ 


= ^^u^yyn,^ + N„ . Any time-domain channel taps 


u-l 


whose magnitude is below the threshold value 7, may be discarded by the 
invention in order to determine the dominant time-taps (time-taps values with 
power above the threshold value) 720. In one embodiment of the invention. 
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the dominant time-taps may be the P time-taps labeled £^ through (p, and may 
represent £ for which \h„,{ef > t . With the dominant time-taps determined, the 

least squares channel estimation method 700 may re-estimate the channel for 
at least one user, only at the dominant time-tap (tap) values 730. For the 
previously described embodiment, the dominant tap values can be.^i through 
£p, whose gains (channel estimations) may be found using a LS procedure 
similar to equation (7), thus terminating 740 the least squares channel 
estimation method 700 

The above-described methods and implementation of adaptive channel 
estimation are example methods and implementations. These methods and 
implementations illustrate but one possible approach for adaptive channel 
estimation. The actual implementation may vary from the method discussed. 
Moreover, various other improvements and modifications to this invention may 
occur to those skilled in the art, and those improvements and modifications 
will fall within the scope of this invention as set forth below. 

The present invention may be embodied in other specific forms without 
departing from its spirit or essential characteristics. The described 
embodiments are to be considered in all respects only as illustrative and not 
restrictive. 
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